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Abstract 

In numerical computations of Einstein's equations for black hole spacetimes, it will be necessary to 
use approximate boundary conditions at a finite distance from the holes. We point out here that 
"tails," the inverse power-law decrease of late-time fields, cannot be expected for such computations. 
We present computational demonstrations and discussions of features of late-time behavior in an 
evolution with a boundary condition. 
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I. INTRODUCTION 



There is at present great interest in the computation of the gravitational waves from the inspiral and merger of 
a pair of mutually orbiting black holes P, 0, Q- To do such computations a solution of the initial value equations 
of general relativity is chosen on some initial spatial hypersurface, and the remaining Einstein equations are used to 
find the spacetime to the future of that initial surface. In principle, one can compute the evolved spacetime only 
in the domain of dependence of that initial surface, which means that the initial surface must be large, many times 
the radius of the initial binary orbit, if the spacetime is to be evolved for several orbital times. The computational 
demands for such a procedure make it unfeasible for the foreseeable future, although pseudospectral codes can help 
considerably in extending the size of the initial hypersurfaceQ. The alternative to a large initial hypersurface is a 
timclike boundary on the computation, typically at some large radius, at which appropriate approximate boundary 
conditions are specified. These boundary conditions are chosen to represent (approximately) the condition that no 
information moves inward through the boundary. One of the problems that workers in this field have recently turned to 
is that of appropriate boundary conditions, especially in connection with the preservation of gauge constraints00,0]. 

One of the features found in evolutions of perturbations in black holes spacetimes is the final latest-time behavior 
of the perturbation fields, a fall off in time t as l/t™ at a constant distance from the hole[^, This was first 
demonstrated for Schwarzschild holes in which case n = 2^ + 3 for a multipole of index £ with compact initial support. 
For Kerr holes such tails also represent the latest time behavior, though there remains some controversy about the 
value of n\Wi Ull Il2|. Most work on these tails has been done within linearized perturbation theory, though some 
computations with self-gravitating spherically symmetric scalar fields have also been carried out [l3lll4Lll5| . 

The purpose of the present paper is to point out that numerical computations with approximate boundary conditions 
cannot be expected to produce the correct late-time tails. A rough intuitive reason for this is that the late-time tails 
are produced by the scattering of radiation due to the curvature of spacetime. This scattering takes place far from 
the hole, and depends on the asymptotic large-distance nature of spacetime curvature. A boundary condition on a 
timelike surface at finite radius means that the asymptotic nature of the distant spacetime does not enter into the 
computation, so that the correct late-time tails cannot develop. 

These ideas can be checked accurately in the Schwarzschild background. In this case linearized perturbations 
(scalar, electromagnetic, or gravitational) can be analyzed into multipolcs ip = ^2 ipem{t, r)Yi m (9, (f>), where t, r, 9 
are the standard Schwarzschild coordinates. The evolution of each multipole is described by a simple wave equation 

dfy(t, r) - d 2 „yj(t, r) + V(r)^(t, r) = . (1) 

Here c = G = 1 and, for simplicity, we have dropped the multipole indices on ipi m (t,r) and on Vi{r). The Regge- 
Wheeler[3 

"tortoise coordinate" r* is defined by r* = r + 2M\n(r/2M — 1), where M is the mass of the Schwarzschild 
hole. In the limit M — > 0, that is, in flat spacetime, r — r* and the potential takes the form 

V Rat = £(£+l)/r* 2 . (2) 

For this potential, the solution to Eq. has a simple familiar form in terms of spherical Bessel functions, and has 
no tails. 

If M ^ there are important differences. For a scalar perturbation, the potential, 
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is typical of that for all perturbation fields. For r* M this potential has the form 

£(£+!) 
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It is the extra In r* /r* term that is in Eq. Q), but missing in Eq. that produces the l/t 2i+3 taiLs[l7j. To check 
this we define a toy potential 
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with an adjustable parameter, A, that allows us to control the size of the extra term. Using this potential we have 
numcr ically evolved Eq. Q with £ = lfis|. The initial data, at time t = 0, for this evolution was a time-symmetric 
Gaussian pulse ip = exp(— 0.1(r*/2M + 4) 2 ). The results, for ip at r* = 400M, are shown as a function of time t in 
Fig. ^ The figure shows the straight lines in the log-log plot that indicate a power-law fall off in time. For all values 
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FIG. 1: Late-time evolved solutions for different values of A, the coefficient of the scattering term in the toy potential described 
in the text. The absolute value of the scalar field ip is shown as a function of time ta.tr*— 400M. Curves are labeled with 
the value of A. 

of A the slope is -5, consistent with the 2£ + 3 rule. (At very late times each of the power-law tails disappears into 
round-off noise.) The results clearly show the sensitivity of the tail to the size of A. For A = 0, of course, the toy 
potential becomes the flat spacetime potential centrifugal, and there is no tail. (The magnitude of the tail increases 
as A is made larger, but this increase slows, and the magnitude of the tail appears to reach a limit.) 

The argument against tails with finite boundary conditions, then, is that they depend on the asymptotic form of 
the potential. For the obvious Sommcrfcld outgoing boundary condition 

ip,t + i>,r* = , at r* = r£ c , (6) 

the interior solution cannot "know" what the potential is for r* > rj^ c , and hence the solution cannot develop a tail 
except at spacetime points whose domain of dependence lies within the outer boundary. 

This argument, of course, is only suggestive. One can rebut it with the claim that a boundary condition could be 
made sufficiently precise so that it encodes the asymptotic form of the potential. A trivial example is the case V = 
and r — r* , for which the general solution, in terms of an arbitrary outgoing function /, and an arbitrary ingoing 
function g, is 



^ = f{t-r)+g{t + r) 



(7) 
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In this case the outgoing boundary condition ipj — ip r — is exact; it constrains g to be zero, the same condition as 
if the boundary were infinitely far away. A less trivial, but less useful case is the M = 0, £ = 1 scalar equation . The 
general solution in this case is 

^ = f'(t-r) + f -^l + 9 '(t + r)- 9 -^, (8) 
r r 

where / describes the outgoing part and g the ingoing part. An exact boundary condition is 

(d r + d t f[r 2 (d r + dt)ij] =0, (9) 

which constrains g (more precisely g ) to vanish. 

The cases in Eqs. J7J) and (JSJ are special. In these cases the radiative part of the solution (the parts without 1/r 
factors) and the nonradiative parts are simply related. The nonradiative part is missing for Eq. JJJ), and for Eq. JHJ the 
nonradiative part is a simple time intergral over the radiative part. This simplicity is related to Huygens' principle. 
Note in particular that in Eq. @ or JSJ if g = and if / has compact support, then the solution will be nonzero only 
for a finite time, and hence cannot have a power-law tail as a generic feature. The nonsimplc non-Huygens solution 
for problems in which tails are generic will not, therefore, have an exact boundary condition that can be expressed as 
a finite number of differentiations as in Eq. @. 



II. EVOLUTIONS WITH APPROXIMATE BOUNDARY CONDITIONS 



We now directly investigate the late-time behavior of solutions with approximate outgoing boundary condition of 
Eq. jnj. The monopole £ — case is studied so that the predicted slowly- falling 1/t 3 tails stay well above the roundoff 
noise. The starting field is a purely outgoing pulse defined on the ingoing null ray t + r* — 0. On this ray, in terms of 

retarded time u = t — r* , the form of tp is specified to be ip — [u(u — 8M)/16Af 2 ] 8 , for < u < 8M, and ip = for 
u < or u > 8M. The initial data at t = then is approximately an outgoing pulse (satisfying ip ir * = —ip,t) confined 
between r* — — 8AT and r* = 0, centered at r* = —AM with a peak value ip = 1 there. The left boundary for all 
I = computations is at r* = — 500M, where the sommerfeld condition -0 ir * = ipj is imposed. The right boundary is 
placed at different locations r£ c , and the condition tp^ r * — —ipj is imposed. 
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FIG. 2: The scalar monopole ip evolving from an initially outgoing pulse near r* = 0. The field ip is plotted as a function of 
time at r* = for two different values of r£ c , at 500M and 200M. The two curves are identical up to time t — 400M. A curve 
proportional to 1/t 3 is included for comparison. 

Figure |21 shows the time profile of ip developing from the initial pulse. The value of ip at the fixed position r* = is 
shown as a function of time for outgoing boundary locations r^ c — 500M (solid curve) and at r^ c = 200M (dot-dashed 
curve). Since the results are shown only up to t = 1000M, the field at r* = 0, shown by the solid curve, has not yet 
been influenced by the interaction of the boundaries at ±500 M with the evolution of the initial pulse from r* = 0. The 
solid curve then shows the "boundaryless" behavior of the field evolving completely within the domain of dependence 
of the initial data. That field goes through quasinormal ringing up to around t = 150M, then around t — 300M 
becomes a power-law tail. Since the plot is a semilog plot, not a log-log plot, the tail is not a straight line, and a 
constant/t 3 curve is provided to illustrate the tail nature of the solid curve. The dot-dashed curve for the r^ c = 200A/ 
outer boundary is identical to the solid curve up to t — 400Af . In particular, from t rj 300Af to t = 400M the field 
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starts to take the form of a power-law tail, but at t = 400Af waves reflected from the r£ c = 200Af boundary arrive at 
r* = and a new oscillatory behavior ensues due to the influence of the boundary. 

A complementary viewpoint on the boundary influence is given in Fig. |3| in which spatial profiles of evolved fields 
are shown at the moment of time t — 500A/. As in Fig. |2 the solid curve corresponds to an outer boundary at 
r* = 500AT, and the dot-dashed curve to an outer boundary at r* = 200A7. The solid curve indicates that at 
t = 500M the unit-height pulse has reached the outer boundary. In the absence of scattering this would be the 
only feature of the plot, but due to scattering there are quasinormal bumps created between r* « 370M and 500M. 
Waves scattered inward also create quasinormal bumps between r* w — 300M and — 500M. From r* w — 300M to 
w 370M the solid curve shows the spatial profile of the tail behavior. The dot-dashed curve shows that reflections of 
the inital pulse off the r* = 200M boundary have reached r* = — 100M and have "contaminated" the spatial profile 
from r* = — 100 A/ to 200Af. The spatial profile from r* = — 100M to 200M shows spatial oscillations similar to the 
temporal oscillations that appear in Fig. [21 

III. BOUNDARY-INDUCED QUASINORMAL RINGING 
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FIG. 3: Spatial profile of the scalar monopole ip evolving from an initially outgoing pulse. Here ip is plotted at t = 500A4 for 
values of r£ c at 500M and 200M. The curves are identical for r* < — 100M where the field has not yet been influenced by 
reflections of the initial pulse off the 200M boundary. 

The nature of the boundary-induced oscillations is made clearer in Fig. which shows the time profiles for outer 
boundaries at lOOAf , 50M, 25Af. (For larger values of r^ c similar oscillations develop at later times.) The late-time 
results for the 50AT and 25Af profiles indicates a constant-period oscillation, and the straight-line envelope of the 
oscillations in the semilog plot indicates an exponential damping. The natural interpretation is that these damped 
oscillations are a new form of quasinormal ringing. The "old" form is the familiar quasinormal ringing of a black 
hole [19l l2fl |21|. like that in Fig. |2] for t/M less than ~ 150. This is a real physical phenomenon associated with the 
black hole spactime. In Fig. 0] we see oscillations that are not physical in this sense, but are numerical artifacts, the 
strongly damped oscillations of a leaky cavity created by the outer boundary and the curvature potential V near its 
peak at r* =0. Some evidence for this is the dependence of the period of the oscillations on r£ c . The longer cavity 
of the 50Af-boundary case creates longer oscillation periods than the cavity of the 25Af -boundary. The less-clcarly 
defined oscillations of the 100M have yet longer periods. 

The leaky-cavity interpretation is supported by the plot in Fig. [S] showing the period T for a full oscillation (the 
width of a pair of "bumps" in Fig.QJ of the late time features, as a function of rj^ c , the location of the outer boundary. 
The uncertainty in the period is around 2%, both due to truncation error in the evolution of ip, and due to the 
extraction of T from the late time results. For r^ c larger than ~ 100M the period quite accurately follows the linear 
relation T/M = —30 + 2.9(r^ c /M ), confirming the view that these late time oscillations are resonances of a cavity 
created by the outer boundary. For r£ c less than ~ 100 M it is not surprising that the details near the peak of the 
curvature potential V would complicate the relationship; for large values of r^ c the details of the curvature potential 
become unimportant. What is at first surprising in the linear relationship is the coefficient 2.9. This means that the 
number of half-wavelengths inside the "cavity" is not an integer; one would naively expect the period T to be equal 
to the cavity length, or an integer multiple of half the cavity length. But our naive expectations are based on the 
simple boundary condition that the field or its normal derivative vanish. In our artificial cavity the outgoing boundary 
conditions does neither. For a constant-frequency oscillation the outgoing boundary condition, in effect, imposes a 
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FIG. 4: Evolution of the t[> monopole for Sommerfeld outgoing conditions on close boundaries. The time profiles are for an 
initial outgoing pulse near r* = 0. Results are shown for outer boundaries at locations r* = 100A/, 50M, and 25M. 
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FIG. 5: The period T for a full oscillation (two "bumps") as a function of the location r^ c of the outer boundary. The straight 
fine T/M = —30 + 2.9(r£ c /M) is the least-squares fit to the results for r^ c = 100M and larger. 



relationship between the field and its normal derivative. Computations with simple toy models have confirmed that 
with such boundary conditions a resonance of the cavity will not contain an integer number of half wavelengths. 

The leaky-cavity viewpoint on the outer boundary suggests a way to quantify the effectiveness of the outer boundary 
condition. The more effective the outer boundary condition is, the more quickly the cavity modes should die out. If 
we characterize the effective reflectivity of the outer boundary as R, and let N represent the number of reflections 
from the outer boundary, then the amplitude should fall off as R N . We can take N to be the time t divided by some 
measure of the time for a reflection. For generality we will take this to be KiT, were k\ is some constant, and T 
is the period for an oscillation. The cavity oscillations should then die off as R t /^ KlT \ or as exp(— t/r), where the 
damping time r is — Tki/ bx(R). The outgoing boundary condition is expected to improve as r£ c increases |22j. so 
it is interesting to make the ansatz R = K 2 M/r^ c . This model predicts t/T = K\j [log(r^ c /M) — log(K 2 )]. Figure|H| 
shows t/T for a range of outer boundary locations r^ c ; the uncertainty here, as in Fig. is no worse than 2% at any 
r^ c . The large- r£ c results in Fig. do show a gradual decrease of t/T with increasing r^ c , as expected. Our heuristic 
model for t is fit to those results by eye, and is plotted in Fig. That fit corresponds to reflectivity 4.8M/r^ c , or 
~ 5% for a boundary at 100M. The plot of the heuristic model gives an appearance of reasonable agreement except 
for small values of r£ c , but with two adjustable parameters in the fit, this agreement can only be said to be weakly 
suggestive. 
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FIG. 6: The time constant r for the envelope of the decaying oscillations, divided by T, the period of the oscillations, as a 
function of the location r£ c of the outer boundary. The thick solid curve shows a heuristic model for r. 

IV. SUMMARY AND DISCUSSION 

We have shown that finite-radius boundary conditions prevent the formation of power-law tails of perturbations in 
black hole spacetimes. In place of a tail, the latest-time feature of a computation will be a new form of quasinormal 
oscillation. Unlike black hole oscillations, these oscillations are not physical phenomena; they arc numerical artifacts 
introduced by the imperfect outgoing condition at the outer boundary. 

In numerical relativity, these oscillation will probably not be a serious practical difficulty. The goal of present 
numerical relativity work is a better understanding of strong field nonlinear dynamics. Numerical codes in the 
forseeable future will not be able to run long enough for the boundary-induced oscillations to appear, nor are they 
likely to be accurate enough to deal with such a weak field phenomenon. 

These limitations of running time and accuracy do not apply when the Lazarus |2^| method is used in a problem 
involving the formation of a final black hole. That method uses the solution computed by a fully nonlinear numerical 
evolution code as initial data for further evolution by black hole perturbation theory. The nonlinear numerical evolu- 
tion would, of necessity, use timelike boundaries, but the numerical cavity oscillations would not develop in the limited 
time for which the evolutions run. In principle, the subsequent Lazarus evolution could be used without boundary 
conditions, i.e., with evolution only within the domain of dependence of the initial data inherited from the fully 
nonlinear numerical code. Such boundaryless Lazarus evolutions would exhibit power-law tails, but the tails would be 
strongly affected by boundary effects contained in the initial data inherited from the nonlinear evolution. In practice, 
Lazarus evolutions are not boundaryless. Rather, to reduce memory requirements, timelike boundary conditions are 
used in the Lazarus perturbation evolutions. These evolutions should contain boundary-induced artifacts of the type 
we have discussed above. But these artifacts would be miniscule, and of no concern for most applications of black 
hole evolution, with or without Lazarus. 
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